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Abstract 

We study the linear and nonlinear stability of relative equilibria in the planar n-vortex 
problem, adapting the approach of Moeckel from the corresponding problem in celestial me- 
chanics. After establishing some general theory, a topological approach is taken to show that 
for the case of positive circulations, a relative equilibrium is linearly stable if and only if it 
is a nondegenerate minimum of the Hamiltonian restricted to a level surface of the angular 
impulse (moment of inertia). Using a criterion of Dirichlet's, this implies that any linearly 
stable relative equilibrium with positive vorticities is also nonlinearly stable. Two symmetric 
families, the rhombus and the isosceles trapezoid, are analyzed in detail, with stable solutions 
found in each case. 
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1 Introduction 

In 2001, Kossin and Schubart J2U] conducted numerical experiments describing the evolution of thin 
annular rings with large vorticity as a model for the behavior seen in the eyewall of intensifying 
hurricanes. In a conservative, idealized setting, they find examples of "vortex crystals," formations 
of mesovortices that rigidly rotate as a solid body. One particular formation of four vortices, 
situated very close to a rhombus configuration, is observed to last for the final 18 hours of a 24-hour 
simulation (see Figure 4(a) of [20J ) . Rigidly rotating polygonal configurations have also been found 
in the eyewalls of hurricanes in weather research and forecasting models from the Hurricane Group 
at the National Center for Atmospheric Research (see [11] or the website [10] for some revealing 
simulations). 

It is natural to explore these rigidly rotating configurations in a dynamical systems setting by 
studying relative equilibria of the planar n-vortex problem. Introduced by Helmholtz [TB] and later 
given a Hamiltonian formulation by Kirchhoff |19| . the n-vortex problem is a widely used model for 
providing finite-dimensional approximations to vorticity evolution in fluid dynamics (3] |29] . The 
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general goal is to track the motion of the vortices as points rather than focus on their internal 
structure and deformation, a concept analogous to the use of "point masses" in celestial mechanics. 
A relative equilibrium is a special configuration of vortices that rigidly rotates about its center of 
vorticity. In a rotating frame these solutions are fixed. Given their persistence in many models of 
hurricanes, analyzing the stability of relative equilibria in the planar n-vortex problem is of great 
interest and may have some practical significance. 

Several researchers have studied the stability of relative equilibria, beginning with the work of 
Lord Kelvin [TBI [T7] . Grobli [T3] and Synge [37] on the well-known equilateral triangle solution 
(the analogue of Lagrange's solution in the three-body problem). If Tj represents the vorticity or 
circulation of the z-th vortex, then placing the vortices at the vertices of an equilateral triangle 
results in a relative equilibrium for any choice of circulations (positive or negative). This periodic 
solution is linearly stable if and only if the total vortex angular momentum L, defined as 

L = J^V^, 

is positive [37]. As we will show, the quantity L turns out to be very important in the analysis of 
stability. We note that both Synge [37] and Aref [4] show, using a special coordinate system, that 
the equilateral triangle solution is actually nonlinearly stable whenever it is linearly stable. 

Another interesting stability result concerns the regular n-gon (n > 4), a relative equilibrium 
only if all vortices have the same strength. It is linearly stable for n < 6, degenerate at n — 7, 
and unstable for n > 8 [38"l [T5"l [2J • By calculating the higher order terms in the normal form of 
the Hamiltonian, Cabral and Schmidt [9l [35] show that in the special case n = 7, the "Thomson 
heptagon" is locally Lyapunov stable. It is possible to extend the regular n-gon relative equilibrium 
by adding an additional vortex of arbitrary circulation v at the center of the ring. Cabral and 
Schmidt show that, for all n > 3, this 1 + n-gon relative equilibrium is locally Lyapunov stable as 
long as v belongs to a bounded interval of values that depends on n [Hj. As was the case with the 
equilateral triangle, for each n, the interval where the 1 + n-gon is linearly stable coincides with 
the one where the configuration is nonlinearly stable. This same phenomenon where linear stability 
and nonlinear stability coincide also occurs in many cases for point vortices on the sphere El |2T] . 
Finally, we mention the recent work of Barry, Hall and Wayne considering the stability of relative 
equilibria containing one "dominant" vortex (z/ >> 0) and n equal strength (r^ = e) but small 
vortices [B]. Following the analogous work of Moeckel in the n-body problem [27], they show that 
a relative equilibrium with e > (e < 0, resp.) in this setup is linearly stable provided the n small 
vortices are situated at a nondegenerate local minimum (maximum, resp.) of a special potential 
function. 

In this work, we follow Moeckel's (2H] approach for studying linear stability of some symmetric 
configurations of relative equilibria in the n-body problem. Working in a rotating frame, the basic 
idea is to find invariant subspaces of the Hessian of the Newtonian potential (multiplied on the left 
by the inverse of a diagonal mass matrix), and use them to determine the linear stability of the 
relative equilibrium. While this is somewhat complicated in the n-body case, it is straight-forward 
for the n-vortex problem, in part because the number of degrees of freedom is just n, rather than 2n. 
Much of our analysis builds from a simple algebraic property of the logarithmic Hamiltonian of the 
planar n-vortex problem: the Hessian anti-commutes with a key rotation matrix. 



2 



One of the interesting conjectures in the n-body problem, due to Moeckel and listed as Problem 
16 in [lj, states that a linearly stable relative equilibrium must be a nondegenerate minimum of the 
Newtonian potential restricted to a level surface of the moment of inertia. For the case of n vortices 
with the same sign, we prove that linear stability is in fact equivalent to being a nondegenerate 
minimum of the Hamiltonian H restricted to a level surface of the angular impulse (the analog of 
the moment of inertia). Since the angular impulse is also a conserved quantity for the n- vortex 
problem, it follows that, in the case of same signed circulations, linear stability implies nonlinear 
stability. In other words, solutions that start near a linearly stable relative equilibrium stay close to 
it for all forward and backward time. This helps explain many of the results described above where 
linear stability and nonlinear stability coincide. We also show that the collinear relative equilibria 
with same signed circulations are all unstable, and that for generic choices of positive circulations, 
a linearly stable relative equilibrium always exists. While previous works have focused on specific, 
often symmetric relative equilibria, our results are more general. 

The paper is organized as follows. In Section 2 we describe the n-vortex problem in rotating 
coordinates, highlight the key properties of the Hamiltonian H, and establish some general theory 
for determining the linear stability of a relative equilibrium. The key quantity L is discussed 
and it is shown that relative equilibria with L = are always degenerate. Section 3 focuses on a 
topological approach to the problem and proves our main result connecting linear stability and being 
a minimum of H restricted to a level surface of the angular impulse. The relation between linear 
stability and nonlinear stability is carefully explained and proven. Finally, in Section 4 we apply 
some of our theory to analyze the linear stability of two families of symmetric relative equilibria 
in the four-vortex problem: a rhombus and an isosceles trapezoid. In both cases, stable examples 
exists for positive circulations. The rhombus case is particularly interesting as it contains two one- 
parameter families of solutions, one of which goes through a pitchfork bifurcation at a parameter 
value expressed as the root of a certain cubic polynomial. Throughout this work, calculations and 
computations were confirmed symbolically using Maple [22] and numerically using Matlab [23J. 



2 Relative Equilibria 



A system of n planar point vortices with vortex strength Ti ^ and positions z\ G 1R 2 evolves 
according to 



Y lZl = JViH = JJ2 



rr 



.-, y Zj — Zi), 1 < i < n 



where 



H = -^TiTjHrij) 



i 1 1 %i Zj | 



i<3 



is the Hamiltonian, J is the standard 2x2 symplectic matrix J 



and V,- denotes the 



1 
-1 

two-dimensional partial gradient with respect to Zj. We will assume throughout that the total 
circulation T = Yli^i * s nonzero. The center of vorticity is then well-defined as c = ^ Yli^i z i- 
Without loss of generality, we take c = 0. Another important quantity is the total vortex angular 
momentum L, defined earlier as L = Ylii<j ^i^j- If a H the vortex strengths are positive, then L > 0. 
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However, when vorticities of different signs are chosen, L < is possible. The case L = turns out 
to be of considerable importance. 

A configuration of vortices which rotates rigidly about its center of vorticity is known as a 
relative equilibrium. 

Definition 2.1. A relative equilibrium is a solution of the form 

Zi (t) = e-^z^O), l<i<n, 
that is, a uniform rotation with angular velocity u^O around the origin. 
The initial positions Zj(0) of a relative equilibrium must satisfy 

-wr jZi (0) = V l H = for each zG{l,...,n}. (1) 

Let z = (zi, z 2 , . . . , z n ) G IR 2n describe the vector of positions. We will typically refer to the initial 
positions zq = z(0) G M 2ra as a relative equilibrium, rather than the full periodic orbit. 

The important quantity / = rj||^|| 2 is the angular impulse, analogous to the moment 

of inertia in the n-body problem. It is well known that / is an integral of motion for the n-vortex 
problem . The system of equations defined by ([I]) can be written more compactly as 

VH(z) +coVI(z) = 0. (2) 

Thus, relative equilibria can be viewed as critical points of the Hamiltonian restricted to a level 
surface of J, with u acting as a Lagrange multiplier. This provides a very useful topological approach 
to studying the problem. 

Note that if z G R 2n is a relative equilibrium, so is rz$ for any scalar r. In this case, the 
angular velocity u must be scaled by a factor of 1/r 2 . Moreover, if e Zq denotes (e J6 zi, . . . , e J9 z n ), 
then e je z is also a relative equilibrium for any 6, with the same angular velocity. Thus, relative 
equilibria are not isolated. When counting the number of solutions, it is customary to fix a scaling 
(e.g., 1 = 1) and identify relative equilibria that are identical under rotation. In this fashion, one 
counts the number of equivalence classes of relative equilibria. We also note that the stability type 
of a particular relative equilibrium is unchanged by scaling or rotation. In particular, if zq is linearly 
stable, then so is any other relative equilibrium equivalent to z Q under scaling or rotation. 

To analyze the stability of relative equilibria, we change to rotating coordinates. Denote M = 
diag{Ti, Ti, . . . , T n , T n } as the 2n x 2n matrix of circulations. Also, let K be the 2n x 2n block 
diagonal matrix containing J on the diagonal. The equations of motion are thus written simply as 
Mz = KVH(z). Under the transformation Z\ h-> e~ wJt Zi Vi, the system is transformed into 

Mz = K(VH(z) +ooMz), (3) 

a uniformly rotating frame with period 2-k/uj. As expected, a rest point of (|3]) is a solution to 
equation ([2]), a relative equilibrium. 
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Remark. Note that the system of differential equations in (J3J) is not written in the usual canonical 
form since the conjugate variables Xj, yi appear consecutively. One could permute the variables and 
use the canonical 2n x 2n matrix J as opposed to K. However, we choose this set-up since it is 
similar to the Newtonian n-body problem and because a key property of the Hamiltonian is more 
readily apparent in these coordinates. 

The • in the following lemma represents the standard Euclidean dot product and D 2 H(z) is the 
Hessian matrix of H. 

Lemma 2.2. The Hamiltonian H has the following three properties: 
(i) VH(z) -z = -L, 

(ii) VH(z) ■ (Kz) = 0, 

(iii) D 2 H(z)K = -KD 2 H(z). 

Proof: Differentiating the identity H(rz) = H(z) — Lln|r| with respect to the scalar r and 
evaluating at r = 1 verifies property (i). If e Jt z denotes (e Jt zi, . . . , e Jt z n ), then the SO (2) symmetry 
of H gives H(e Jt z) = H(z) Vt. Differentiating this identity with respect to t and evaluating at 
t = shows property (ii). To prove the third property, notice that differentiation of (i) implies 
-KD 2 H(z)z = KVH(z), while differentiation of (ii) yields D 2 H(z)K z = KVH(z). Since z is 
arbitrary, these two equations prove (iii). □ 



Using property (i) and the fact that / is homogeneous of degree 2, equation ^ implies that the 
angular velocity of a relative equilibrium is u> = L/(2I). If all the vortex strengths are positive, then 
u > is assured. If L = 0, then we must also have I = in order to have a relative equilibrium. 
Since dl/dt = —2VH(z) ■ (Kz), property (ii) verifies that the moment of inertia is a conserved 
quantity in the planar n-vortex problem. 

The fact that the Hessian of the Hamiltonian and the rotation matrix K anti-commute will 
be profoundly important in the stability analysis of relative equilibria. For positive vorticities, it 
implies a complete factorization of the characteristic polynomial into even quadratic factors. We 
note that property (iii) holds true for any degree zero homogeneous potential with SO (2) symmetry. 

2.1 Linear Stability 

We will assume throughout that zq G IR 2n represents a particular relative equilibrium. Linearizing 
system ^ about z gives the stability matrix 

B = K(M- 1 D 2 H(z ) +uT), 

where I is the 2n x 2n identity matrix. Here we have used the fact that K and M _1 commute. 
Through a permutation of the variables, our stability matrix B is similar to the matrix used in [U] . 
The eigenvalues of B determine the linear stability of the corresponding periodic solution. Since 
the system is Hamiltonian, they come in pairs ±A. To have stability, the eigenvalues must lie on 
the imaginary axis. The eigenvalues of B can be described simply in terms of the eigenvalues of 
M~ 1 D 2 H(zq). The next two lemmas make this statement precise. 
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Lemma 2.3. The characteristic polynomials of D 2 H(z$) and M^ 1 D 2 H(zq) are even. Moreover, 
for each matrix, v is an eigenvector with eigenvalue \x if and only if Kv is an eigenvector with 
eigenvalue —fi. 

Proof: We give the proof for M~ l D 2 H(zo); the proof for D 2 H(zo) is basically identical. Note that 
det(K) = 1 and K 2 = —I. Let q(X) = det(M~ 1 D 2 H(z ) — XI) be the characteristic polynomial for 



M 1 D 2 H(zq). Then, using property (iii) of Lemma 2.2, we have 



q(X) = det( y K)det(M- 1 D 2 H(z ) - XI) det(K) 
= det( K(M~ 1 D 2 H(z )K + XI) 
= det( M^ 1 D 2 H(z ) + XI) 

= g(-A). 

We also have the following sequence of implications: 

M- 1 D 2 H(z )v = fiv <^> KM- 1 D 2 H(z )v = fiKv 

<^ -M- l D 2 H(z )Kv = fiKv 
<^=> M- 1 D 2 H(z )(Kv) = -fi(Kv), 

which verifies the second part of the lemma. □ 



Lemma 2.4. Let p(X) denote the characteristic polynomial of the stability matrix B. 

(a) Suppose that v is a real eigenvector of M _1 D 2 iJ(z ) with eigenvalue \i. Then {v,Kv} is a 
real invariant subspace of B and the restriction of B to {v, Kv} is 

fi-uA () 
fi + u J ' y } 

Consequently, p(X) has a quadratic factor of the form X 2 + u 2 — /i 2 . 

(b) Suppose that v — V\ +1V2 is a complex eigenvector of M~ 1 D 2 H(zq) with complex eigenvalue 
fi = a + i (3. Then {v 1, v 2 , Kv\, Kv 2 } is a real invariant subspace of B and the restriction of 
B to this space is 

a-u (3 

—/3 a — u , > 

a + co f3 W 

-p a + u 

Consequently, p(X) has a quartic factor of the form (A 2 + u 2 — tt 2 )(X 2 + uj 2 — Jl 2 ). 



Proof: First, suppose that M 1 D 2 H(z )i 
Bv = (/x + u)Kv and B(Kv) — (// — u)v 



[iv, with v real. Using Lemma 2.3, this implies that 



matrix Q is A 2 + u 2 - jj, 



verifying matrix Q. The characteristic polynomial of 
2 and consequently, this quadratic must be a factor of p(X). 
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Second, suppose that M 1 D 2 H{zq)v = fav, with v = v±+i v 2 and \i = a+i (3. Using property (iii) 
of Lemma 2.2, this implies that Bv\ = (a + u)Kv\ — (3Kv 2 , Bv 2 = (a + iS)Kv 2 + 0Kvx, B(Kvx) 



(a 



uj)vi — (3v 2 and B(Kv 2 ) 
polynomial of matrix ^ is (A 2 



— tu)v 2 + (3vi, which verifies matrix ([5]). The characteristic 



fjL ) and consequently, this quartic, which has 



real coefficients after expanding, must be a factor of p(X). 



□ 



Lemma [2T4] suggests a simple strategy for determining linear stability of a relative equilibrium zq. 
If enough eigenvectors of M~ 1 D 2 H(z ) can be found, then the corresponding eigenvalues fa deter- 
mine stability. If fa G K, then \fa\ < \u\ is required for linear stability. If fa 6 C, then Re(fa) = 
is required for linear stability. 

Write the relative equilibrium as z = (z±, z 2 , ■ ■ ■ , z n ) e M? n and let Zij = {z{ — Zj)/rij. Direct 
computation shows that 

'An A 12 ■ 



D 2 H(z ) 



where Aj,- is the 2x2 symmetric matrix 



A n \ A. 



ii2 



A ln 

A 



(6) 



A 



L ' rj [/-2 % 4] ifz^j, 



r 



An 



(7) 



Note that Ay = A,j and that for % ^ j, 



A 



r r 



(Vi - Vj) 2 - fa - x j) 2 ~2(xi ~ Xj)(yi - y s ) 
-2(xi - Xj)(yi - Vj) (xi - xj) 2 - {yi - yj f 



The fact that J commutes with each Ay gives another proof that D 2 H(zq) and K anti-commute. 

Any relative equilibrium in the planar n-vortex problem will always have the four eigenvalues 
0, 0, We will refer to these eigenvalues as trivial. They arise from the symmetry and integrals 
of the problem. Due to conservation of the center of vorticity, the vectors s = [1, 0, 1, 0, ... , 1, 0] 
and Ks are in the kernel of D 2 H(zq). This is easily verifiable from the structure of D 2 H(zq). 
Consequently, \x = in Lemma 2.4| and B has the eigenvalues ±ui. In addition, due to the fact that 
relative equilibria are not isolated under rotation, the vector Kz is in the kernel of B. This can be 
verified analytically by differentiating property (i) of Lemma 2.2 and substituting in equation ([2]). 
This gives the important identity 



M^D^Hizo) z 



uz , 



which holds for any relative equilibrium z . Equation (|8]) together with property (iii) of Lemma 2.2 
then shows that BKzq = 0, as expected. 

For a given relative equilibrium zq, let V = spanj^o, Kzq}. This is an invariant subspace for B 
and the restriction of B to V is given by 

_ 
2cu 



7 



Thus, a relative equilibrium is always degenerate in our coordinates. The off-diagonal term repre- 
sents the fact that scaling z gives another relative equilibrium but with a scaled angular velocity u. 
We follow Moeckel's approach in [2S] and define linear stability by restricting to a complementary 
subspace of V. 

Recall that M = diag{Ti, F\, . . . , F n , T n }. It is important to note that, unlike the ra-body 
problem, M will be indefinite if two vortices have circulations of opposite sign. We say that the 



vectors v and w are M-orthogonal if v T Mw 
complement of V, that is, 



0. Next, let V 1 - C 



i> In 



denote the M-orthogonal 



{w G R 2n :w 1 Mv = Vf G V}. 



Lemma 2.5. The vector space V 1 - has dimension 2n — 2 and is invariant under B. If L ^ 0, then 
VDV ± = {0}. 



Proof: Note that w G V 1 - if and only if w is in the null space of the 2 x 2n matrix 

zlKM 



(9) 



Since the rank of this matrix is always two, the dimension of the kernel, and hence the dimension 
of V ± , is 2n - 2. 

Next, suppose that w G V 1 - and consider the vector Bw. Using equation ([8j and property 
(iii) of Lemma 2.2, we have that B t Mzq = and B T MKz$ = 2uMzq- This in turn implies 
that (Bw) T Mz = and (Bw) T MKz = 2u(w T Mz ) = 0. Thus, Bw G V L and V L is invariant 
under B. 

Finally, suppose that v G V fl V 1 - . Write v = az + bKz Q for some constants a and b. Since v is 
also in V x , we have 



(az + bKz ) T Mz 
(az + bKz ) T MKz 



0. 
0. 



(10) 
fill 



But z T MKz = for any z G M 2n and since zq is a relative equilibrium, z^Mzq = 21 = L/u. Thus 
) reduces to aL = and 
and v = 0, as desired. 



equation (10) reduces to aL = and equation (11) becomes bL = 0. Since L ^ was assumed, we 
have a = b 



□ 



As long as L ^ 0, Lemma 2.5 allows us to define linear stability with respect to the M-orthogonal 
complement of the subspace V. Instead of working on a reduced phase space (eliminating the 
rotational symmetry), we will stay in the full space and define linear stability by restricting B 
to V ± . 

Definition 2.6. A relative equilibrium z always has the four trivial eigenvalues 0, 0, ±ui. We 
call z nondegenerate if the remaining 2n — 4 eigenvalues are nonzero. A nondegenerate relative 
equilibrium is spectrally stable if the nontrivial eigenvalues are pure imaginary, and linearly stable 
if, in addition, the restriction of the stability matrix B to V x has a block-diagonal Jordan form 

"OA 



with blocks 



-A o 
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2.2 The Special Case L = 



Under this definition, a relative equilibrium (u 7^ 0) with L = will always be degenerate. This 
fact is interesting in its own right; changes in stability and bifurcations are expected to occur when 
L = 0. This is precisely what transpires for the equilateral triangle solution as well as for the 



rhombus B family (see Section 4.2) 



Theorem 2.7. Suppose that zq is a relative equilibrium for a choice of circulations with L 
Then two of the nontrivial eigenvalues for zq are zero and zq is degenerate. 



0. 



Proof: Suppose that L = 0. By Lemma |2.5[ V has dimension 2n — 2 and is invariant under B. 
If Zq satisfies equation ^ for some u 7^ 0, then we must have 21 = ZqMzq = 0. This leads to the 
rather perverse setting where V = span{z , Kzq} is contained in V x , its M-orthogonal complement. 
Using matrix (JoJ) as a guide, consider a basis £ of lR 2n of the form 

£ = {Mz ,MKZq, Z ,Kz Q ,S,Ks,Wi,. . . ,W2n-6} 

where the Wi represent a set of 2n — 6 linearly independent vectors that complete a basis for V x . 
Note that the first four vectors in £ are mutually orthogonal. 

To find the stability matrix B written with respect to £, write 



BMzq = aiMz + a 2 MKz Q + a 3 z + o^Kzq + 
BMKz = biMz + b 2 MKz + b 3 z + b A Kz + ■ 



(12) 
(13) 



o 2 Zq M zq, respectively. Since Zq M z > 0, we 



Taking the standard dot product of both sides of equation (12) with the vectors Mz and MKzq, 
respectively, gives = aiZ F M 2 z and 2uZqM 2 z 
have ai = and a 2 = 2uj. 

the vectors Mzo and MKzo, respectively, gives = b\z^M 2 ZQ and = b 2 ZQM 2 zo, respectively. It 
follows that bi = b 2 = 0. Thus, the stability matrix B written with respect to the basis £ has the 
form 



Similarly, taking the dot product of both sides of equation ( JX3[ ) with 

b 1 z^M 2 z and 





2u 









2u 





-uj 







* 1 u 

and zq has two additional zero eigenvalues in addition to the expected eigenvalues 0, 0, ±coi. □ 



3 A Topological Approach 

In this section we present our main results, focusing specifically on the case Tj > Vj. Recall that 
the angular velocity u = Lj (21), so u > is assured for this case. We also have that M is positive 
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definite. Thus we can define an inner product by 



< v, w > — v Mw. 



Let 8j = [[^]], where [[•]] is the greatest integer function. Denote {e-,} as the standard basis 
vectors of IR 2n and let pj = (1/ wT^TJe^. The vectors {pi,P2, ■ ■ ■ ,P2n} form an M-orthonormal basis 
of IR 2n , that is, < Pi,Pj >= where 5ij is the usual Kronecker delta function. The key matrix 
M~ 1 D 2 H(zq) is symmetric with respect to this basis. Consequently, all its eigenvalues are real 
and it has a full set of M-orthogonal eigenvectors. By Lemma 2.3[ the eigenvalues are of the form 
±//j. It follows that the characteristic polynomial of the stability matrix B factors completely into 
quadratic factors of the form A 2 + uj 2 — /i 2 . 

Theorem 3.1. IfTj > Vj ; then a relative equilibrium zq is linearly stable if and only if \p,j\ < u 
for all eigenvalues fij of M~ 1 D 2 H(zq). 

Proof: For the case Tj > Vj, all the eigenvalues of M~ 1 D 2 H(zq) are real and come in pairs ±/Xj. 



Part (a) of Lemma 2.4 now applies repeatedly and the characteristic polynomial of B factors as 



n-2 



p(A) = A 2 (A 2 + u; 2 )n(A 2 + ^ 2 -/! 



i=i 



If zo is linearly stable, then u 2 — /i 2 > is required for all j. Conversely, if \fij\ < u> for all j, we 
are assured that all the nontrivial eigenvalues are nonzero and lie on the imaginary axis. Thus, z Q 
is nondegenerate and since M~ 1 D 2 H(z ) has a full set of eigenvectors of the form {vj,Kvj}, the 
restriction of the stability matrix B to V x has a block-diagonal Jordan form with blocks 







uj 2 - n 2 



Consequently, z is linearly stable. 



□ 



Remark. 1. For the case of positive circulations, if Zq is linearly stable, then the angular velocity 
of each component of the linearized flow is less than or equal to the angular velocity of Zq. 
If we assume further that fij ^ Vj, then, working on a reduced space that eliminates the 
50(2) symmetry, the Lyapunov center theorem applies to produce a one-parameter family of 
periodic orbits emanating from z 



2. Theorem 3.1 also applies to the case Tj < Vj, although here we have u < so the condition 



for linear stability becomes \fj,j\ < —uj. 

Due to the structure of the Jordan form, it is not possible to have Krein collisions for this 
case. In other words, stability can only be lost by having a pair of pure imaginary eigenvalues 
merge together at and then move onto the real axis; eigenvalue pairs that collide on the 
imaginary axis will pass through each other. With positive circulations, the linear stability 
and spectral stability of z are actually equivalent concepts. 
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Next we follow a topological approach to prove our main result. Recall that I = | Y^l=\ r £ 1 1 ^ 1 1 2 
and that a relative equilibrium is a critical point of H restricted to the level surface I = Iq. Assuming 
all the vortex circulations are positive, I = Iq is an ellipsoid or topological sphere in M 2n . If we 
restrict H to this sphere, it turns out that the local (nondegenerate) minima of H are exactly the 
linearly stable relative equilibria. 

First we recall some important definitions. Suppose that zq is a relative equilibrium with I(zq) = 
Iq. Let S denote the surface I = Iq and let 7 = 7(0) be the closed curve e J9 z of relative equilibria. 
The curve 7 (a circle) lies on S and H is constant on 7. The tangent space of S at z , denoted T Z0 S, 
is simply the M-orthogonal complement of Zq. In other words, T ZQ S — {y £ lR 2ra : v T Mzq = 0}. 
The tangent vector to 7 at z is just Kzq and since < Kzq, zq >= 0, it lies in T ZQ S, as expected. 

Define the Lagrangian function G(z) = H(z) + col(z). By equation (|2]), zq is a critical point 
of G. The type of critical point is determined by examining the quadratic form q Zo (v) associated to 
the Hessian matrix of G [32] • To work with this quadratic form, it is easier to remain in IR 2n rather 
than use local coordinates. We have 

q zo (v) = v T D 2 G(z )v = v T (D 2 H(z Q )+uM)v, veT Z0 S. 

The nullity of q Zo is the dimension of the null space of D 2 G and the index (or Morse index) of q ZQ is 
the dimension of the maximal subspace for which q ZQ (y ) < 0. The connection between the stability 
matrix B and the quadratic form q ZQ is given by the equation —MKB = D 2 G. It follows that z is 
a nondegenerate relative equilibrium if and only if the nullity of q ZQ is one. Therefore, treating zq 
as a critical point of H restricted to S, we will call zq a nondegenerate critical point if the nullity 
of q Zo is one. In this case, since H is constant along 7, if the index of q Zo is zero, then z is a 
nondegenerate minimum of H restricted to S. 

Theorem 3.2. IfTj > Vj ; then a relative equilibrium zq is linearly stable if and only if it is a 
nondegenerate minimum of H subject to the constraint I = Iq. 

Proof: Suppose that z is linearly stable. Since Tj > Vj, the matrix M~ 1 D 2 H(zq) has an 
M-orthonormal set of eigenvectors of the form {vj,Kvj} with eigenvalues ±/ij. Let us choose 
v\ = (1/a/2/ )zq and v% = (1/vT )[1, 0, . . . , 1, 0] T . The vectors {Kvi, t>2, Kv%, . . . , v n , Kv n } form an 
M-orthonormal basis for T ZQ S. Note that 

D 2 G(z )vj = (fij + u)Mvj and D 2 G(z )Kvj = (-fij + u)MKvj. 

Given any vector v £ T Zo S, write it as v — YTj=2 c j v j + Sj=i dj{Kvj). We then compute that 

n 

q zo {v) = u{c 2 2 + d 2 2 ) + ^2[c 2 {fij + u) + d 2 {-fi j 

J=3 



By Theorem 3.1, we have — cu < jij < uo Vj. It follows that q zo (v) > for all v £ T Z0 S. Moreover, 
Qz {v) = if and only if v £ span{i^z }, a direction tangent to the curve of critical points 7. Thus, 
the nullity of q Zo is one and the index is zero. Hence, zq is a nondegenerate minimum of H subject 
to the constraint I = Iq. 

Conversely, suppose that zq is a nondegenerate minimum of H restricted to S and by contra- 



diction, suppose that z is not linearly stable. From Theorem 3.1 there must exist an eigenvector 
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v of M^ 1 D 2 H(z ) with eigenvalue \i such that v $ span{z , Kzq] and /i > u. (If // < 0, then 
replace v with .fTf.) We can also scale v so that v T Mv = 1. We then have Kv G T^S" and 
q Z0 (Kv) = u — < 0. But this contradicts the fact that Zq is a nondegenerate minimum of -ff 
subject to the constraint I = I . □ 



Remark. As mentioned in the introduction, the discovery of this theorem was motivated by the 
related conjecture of Moeckel's in the n-body problem (Problem 16 in [T]). Unlike the case of 
n vortices, the converse direction is clearly false for the n-body problem. For example, Lagrange's 
equilateral triangle solution in the three-body problem is a minimum for any choice of masses but 
is linearly stable only for a small set of mass values. 

For positive circulations, the existence of a minimum of H restricted to I — I® is clear. As 
rij — > 0, H — > oo and since 5* is compact, H restricted to / = Iq is bounded below. It follows that 
H restricted to S has a global minimum, and therefore a relative equilibrium zq exists. However, 
for special choices of the circulations, Zq may be degenerate. For example, the relative equilibrium 
consisting of a regular n-gon (equal strengths) and a central vortex of circulation r is degenerate 
for certain values of To- The number of degenerate cases grows with n (see j9]). Thus, the existence 
of a linearly stable relative equilibrium is assured for most choices of circulations. It turns out that 
this stable relative equilibrium is never a collinear configuration. 

Given n positive circulations, there exist exactly n!/2 collinear relative equilibria, one for each 
ordering of the vortices (the factor of 1/2 occurs because configurations identical under a 180° 
rotation are identified). This follows by generalizing a well-known result of Moulton's for the n- 
body problem [281 EQ]- It can be shown that each collinear relative equilibrium is a minimum of 
H restricted to the sub-manifold of collinear configurations on S. However, in directions normal 
to the collinear manifold, the value of H decreases. In other words, collinear relative equilibria are 
minima restricted to the collinear case but saddles in the full planar problem. Intuitively, as the 
vortices move off the line, the mutual distances increase and the value of H decreases. 

A rigorous proof of this fact is subtle. In [32], Palmore states that the collinear relative equilibria 
have index n — 2 in the planar n- vortex problem, although no proof is given. The same result is true 
for the planar n-body problem and follows from a very clever argument due to Conley, articulated 
by Pacella in [31]. It is straight-forward to check that the proof given in [21] generalizes to the 
planar n-vortex problem. The basic idea is to view the linear flow corresponding to the normal 
directions as a flow on the space of lines through the origin. One can find a "cone" about the trivial 
eigenvector Kzq that is mapped inside itself. This shows that the eigenvalues of M~ 1 D 2 H(zo) 
corresponding to the normal directions are all less than —oj, which means the index of q zo is n — 2. 
The reader is encouraged to see [21] or the lecture notes of Moeckel [23] for more details. Since 
collinear relative equilibria with positive vorticities are never local minima, we have the following 
result. 

Corollary 3.3. Any collinear relative equilibrium with positive vorticities is unstable. All the non- 
trivial eigenvalues come in real pairs ±a_j. 

The study of collinear relative equilibria with circulations of mixed sign is far more subtle, 
as discovered by Aref in his comprehensive study of three- vortex relative equilibria [5J. In this 
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case, it is possible to have linearly stable solutions. We note that Aref's result for the case of 



three collinear vortices with positive circulations is in agreement with Corollary |3.3| Since collinear 
relative equilibria are saddles, there must exist a strictly planar global minimum for H restricted 
to S. This yields another nice corollary of Theorem |3.2 



Corollary 3.4. For a generic choice of positive circulations, there exists a non- collinear, linearly 
stable relative equilibrium. 



Remark. Corollary 3.4 is known to be false for the n-body problem. In fact, if n > 24,306, then 



all of the equal-mass relative equilibria are unstable 



Since both H and / are conserved along solutions of the planar n-vortex problem, Theorem 3J2 
actually has deeper implications involving nonlinear stability. Suppose that z is a linearly stable 
relative equilibrium. Let 7 = e~ ujJt z be the periodic solution corresponding to z . Each point on 
7 satisfies equation ^ and is locally a nondegenerate minimum of H restricted to I — Iq — I(zq). 
Modifying a classical result of Dirichlet's [12J, we will show that solutions starting close to zq will 
stay close to 7 (as a set) for all time (forward and backward). 

To be precise, let 0(t, Co) denote the solution to the differential equation ^ with initial condi- 
tion Co- Define the distance between a point Co an d a compact set Z to be 

d((o,Z) = min{||Co - z\\ : z G Z}. 

We will call the relative equilibrium z nonlinearly stable or simply stable if for any e > 0, there 
exists a 5 > such that ||Co — ^0 1 1 < $ implies d((j)(t, Co), 7) < e for all tel. Since equation ^ is 
the planar n-vortex problem in a rotating frame, solutions starting in a neighborhood of a stable 
relative equilibrium will stay close in the original coordinates as well. 

Theorem 3.5. Suppose Ti > Vz. Then any linearly stable relative equilibrium is also nonlinearly 
stable. 

Proof: We follow the related proof given in [21]. Suppose that z is a linearly stable relative 



equilibrium and define 7 as above. Fix e > 0. Using property (ii) of Lemma 2.2, it is straight- 
forward to check that G = H + col is a conserved quantity for system (|3]). By adding a constant to 
the Hamiltonian H, we can assume, without loss of generality, that G(z ) = 0. Since zq is linearly 



stable, Theorem 3.2 implies that G is strictly positive in a neighborhood of z$ except along 7 where 
it vanishes. The same is true for any point on 7. Let Y v = {z G IR 2n : < d(z,j) < f]} be the open 
tube of radius 77 surrounding 7. Since 7 is compact, there exists an 77 > such that G is positive 
on Y v . 

Let e* = min{e, 77} and let c* = min{G(z) : z G dY e *}. Thus, c* > 0. Since G is continuous away 
from collisions, there exists a 5- neighborhood of Zq where G < c*. For any Co i n this neighborhood, 
we have G(<p(t, Co)) < c* for all t. If there exists some time t* where d(<p(t*, Co), 7) = e*, then 
(j(0(t*,Co)) > c *, which is a contradiction. Therefore, we have d((f>(t, Co), 7) < e* < e for all t, as 
desired. □ 
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4 Some Symmetric Examples 



In this section we apply some of the general theory established above to some specific examples 
of relative equilibria in the planar n-vortex problem. The examples considered are all symmetric, 
which allows for a simpler analysis of the stability. 



4.1 The Equilateral Triangle 



We begin with the well-known equilateral triangle solution in the three-vortex problem. Placing 
three vortices of any strength at the vertices of an equilateral triangle yields a relative equilibrium. 
According to Newton [29] , Lord Kelvin was the first to conclude that the equilateral triangle solution 
was stable for the case of identical circulations [HI [17]. The general case was handled later by 
Synge [37] who used trilinear coordinates to show that it was linearly stable if and only if L > 0. 



We reproduce this result here using our setup and applying Lemma |2.4 

Given any three circulations ri,T 2 
and c 
and z Q 



r 3 , let z\ = (1,0), z 2 = (-l/2,V3/2),f 3 = (-1/2,-^/2) 
F Tlii^i^'i- Then the vortex positions = z \ — c will have center of vorticity at the origin 
Z\,Z2,Zs) is a relative equilibrium with angular velocity tu = T/3. Using equations ^ 



and ([7]), we compute that M 1 D 2 H(z ) is given by 



1 
6 
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-r\ 




ri - 2r 3 
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-(ri - 2r 3 ) 





-2r 3 


-ri 




2r 2 





ri - 2r 2 




-VST, 







-2r 2 




-(r 1 -2r 2 ) 



We know that M~ 1 D 2 H{zq) has the four trivial eigenvalues 0,0, cu, —u. Let ±n be the remaining 
two eigenvalues, where [i may be real or pure imaginary. Instead of searching for an eigenvector, we 
follow Moeckel's approach in [26] and use a formula involving the trace of the square of M^ 1 D 2 H(zq) 
to obtain the remaining eigenvalues. The characteristic polynomial of M~ 1 D 2 H(zq) is 



A 6 + aA 4 + 6A 2 = A 6 -(/i 2 + W 2 )A 4 + /A; 2 A 2 



(14) 



On the other hand, the coefficient a in (14) can be found using a = |[(Tr(C)) 2 — Tr(C 2 )] where 
C = M~ 1 D 2 H(zq) and Tr(*) is the trace of a matrix. This formula applies nicely here because 
Tr(C) = 0. A straight-forward computation gives Tr(C 2 ) 



4p 2 
9 1 



|L. Then we have 



2T 2 L 
"~9~ + 3 



iv 2 + co 2 ) 



E! 

~9~ 



9 



L 
3" 



Expanding fi 2 in terms of the circulations Tj yields 
1 



9 



(r 2 + r 2 + r 2 - r^, - r x r 3 - r 2 r 3 ) 



36 



[(r 1 + r 2 -2r 3 ) 2 + 3(r 1 -r 2 ) 2 ] 
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which shows that /i > with eq ualit y if and only if the three circulations are equal. Note that 
2 



UJ 



[i 2 = L/3. Applying Lemma 2.4, we see that the nontrivial eigenvalues of z Q are 
Therefore, the equilateral triangle is linearly stable when the circulations satisfy L > 0, unstable 
when L < and degenerate when L = 0. This agrees with the result of Synge [37J. We also point 
out that for the case of equal circulations, fi = and the nontrivial eigenvalues are ±ui. In this 
case, the Floquet multipliers for the periodic solution (in the full phase space) are all equal to one. 

4.2 The Rhombus Families 

For the planar four-vortex problem, Hampton, Santoprete and the author showed that there exists 
two families of relative equilibria where the configuration is a rhombus [13]. Set Ti = T 2 = 1 
and T 3 = T 4 = m, where m G (—1,1] is a parameter. Note that L = m 2 + 4m + 1 vanishes 
when m = —2 + a/3 ~ —0.268, so this is an expected bifurcation value. Place the vortices at 
z\ = (1,0), zi = (— 1,0), ^3 = (0, y) and z^ = (0, — y), forming a rhombus with diagonals lying on 
the coordinate axis. This configuration is a relative equilibrium provided that 



if 



X - [p ± v / /3 2 + 4m) , /3 = 3(1 - m), (15) 

or equivalently, 



3|/ 2 - y 4 
3y 2 -l 

The angular velocity is given by 



m = — 2 — r- ( 16 ) 



m 2 + 4m + 1 1 2m 
w = -TTz ^— = - + 



2(1 + my 2 ) 2 y 2 + 1 

where the last equality holds due to the relation between y and m. The cases m < -1 or m > 1 
can be reduced to this setup through a rescaling of the circulations and a relabeling of the vortices. 

There are two solutions depending on the sign choice for y 2 (see Figure [l] for an example). 
Taking + in (15) yields a solution for m e (—1, 1] that always has u > 0. We will call this solution 



rhombus A. Taking — in (15) yields a solution for m e (—1, 0) that has to > for m G (—2 + \/3, 0), 
but uj < for m G (—1, —2 + a/3). We will call this solution rhombus B. At m = —2 + a/3, rhombus 
B actually becomes an equilibrium point (uj = or VH(zo) = 0), while rhombus A is a degenerate 



relative equilibrium (since L = and uj ^ 0, Theorem 2.7 applies). 

It is straight forward to check that for either sign, y 2 is a monotonically decreasing function 
of m. For rhombus A, y 2 has a supremum of 3 + 2a/2 at m = -1 and a minimum of 1 at m = 1 
(the square configuration). For rhombus B, y 2 has a supremum of 3 - 2a/2 at m = —1 and an 
infimum of at m = (vortices 3 and 4 collide). At the bifurcation value m = —2 + a/3, we have 
y 2 = 2 + a/3 for rhombus A and y 2 = 7 — 4a/3 for rhombus B. The graph of y as a function of m 
for each family is shown in Figure [2j 

As explained in [H], the rhombus B family undergoes a pitchfork bifurcation at m = m* ~ 
—0.5951, where m* is the only real root of the cubic 9m 3 + 3m 2 + 7m + 5. Asm increases through 
m*, rhombus B bifurcates into two convex kite configurations. We show here that one pair of 
eigenvalues changes from pure imaginary to real as m passes through m*, as expected at such a 
bifurcation. 
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Figure 1: The two distinct rhombi relative equilibria when m = —0.2. The solutions rotate in 
opposite directions. The configuration on the left (A) is stable; the one on the right (B) is unstable. 



The symmetry of the configuration allows for an exact computation of the eigenvectors and 
eigenvalues of M~ 1 D 2 H(zq). The calculations were done by hand and confirmed using Maple [22J. 
The 2x2 off-diagonal blocks of D 2 H(zq) are given by 
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4 
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2y 



-1 





A 
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m 
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2y 
i-y 2 



A 
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m 







ii 



(y 2 + iy 



y 2 - 1 -2y 
-2y 1 - y 1 



A23 = A14 and A24 = A\^. Due to this nice structure, it is possible to determine the nontrivial 
eigenvectors of M~ 1 D 2 H(zq). They are real for all m and given by 

vi = [my,Q, -my,0, 0, -1,0, 1] T , 
v 2 = [m,0,m,0, -1,0,-1, 0] T , 

Kv\ and Kv%. One can check that these four vectors are M-orthogonal to each other and to the 
trivial eigenvectors. For m > 0, the vector v\ corresponds to a perturbation that keeps the symmetry 
of the configuration but moves vortices 1 and 2 away from the origin while pushing vortices 3 and 4 
toward the origin. This is precisely the kind of perturbation that causes instability for the 1 + n-gon 
relative equilibrium (particularly when n is even) in the Newtonian problem [33] . 
The corresponding eigenvalues for V\ and respectively, are 



/'i 



/'2 



2 (y2 + 1)2 

2(m + 1)(1 -y 2 
(y 2 + I) 2 



(m(y 2 - 1) + 2) 



7y 4 - 18y 2 + 7 



2(y 2 + l)(3y 2 -iy 

2(y 2 -l)(y 2 + 2y-l)(y 2 -2y-l) 
(ya + l)W-l) 
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Figure 2: The value of the position variable y as a function of m. The upper curve corresponds to 
rhombus A while the lower curve represents rhombus B. 

We note that the second expressions for in and /i2 are well-defined since y 2 > 1 > 1/3 for rhombus 
family A and y 2 < 3 — 2y / 2 < 1/3 for rhombus family B. Since we have a full set of M-orthogonal 



eigenvectors for M D H(zq), we can apply Lemma 2.4 to obtain the four nontrivial eigenvalues 



Using equation (16), the key quantities for determining linear stability are 



2 _ 2 ~% 4 - % 2 + l)(3y 4 -2y 2 + 3) 
U Ml ~ (y 2 + l) 2 (3y 2 -l) 2 

^~ A = 4^+1^-^ ^) = 3y 3 -15, 2 + 41,-5. 
Both of these quantities need to be positive to have linear stability. 

Theorem 4.1. 1. Rhombus A is linearly stable for —2 + \/3 < m < 1. At m = -2 + y/3, the 
relative equilibrium is degenerate. For — 1 < m < —2 + \f2>, rhombus A is unstable and the 
nontrivial eigenvalues consist of a real pair and a pure imaginary pair. 

2. Rhombus B is always unstable. One pair of eigenvalues is always real. The other pair of 
eigenvalues is pure imaginary for — 1 < m < m* and real for m* < m < 0, where m* w 
—0.5951 is the only real root of the cubic 9m 3 + 3m 2 + 7m + 5. At m = m* , rhombus B is 
degenerate. 

Proof: Consider u 2 - /i 2 as a function of y. Since 3y 4 - 2y 2 + 3 = 3(y 2 - 1/3) 2 + 8/3 > 0, the 
sign of u 2 — /i 2 is determined by the sign of y A — 4y 2 + 1. Note that this quartic has roots when 
y 2 = 2 ± v3. It follows that the pair of eigenvalues corresponding to {v\, Kv{\ are pure imaginary 
if 2 - < y 2 < 2 + \/3, real if < y 2 < 2 - V3 or y 2 > 2 + \/3, and when y 2 = 2 ± y/3. Since 
y 2 is monotonically decreasing for both families, we can easily translate these statements in terms 
of m. For rhombus A, we have that 1 < y 2 < 2 + \/3 for 1 > m > —2 + \/3, so the eigenvalues 
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are pure imaginary in this interval of m-values. When —2 + \/3 > m > — 1, y 2 > 2 + \/3, and 
the eigenvalues form a real pair. At m = —2 + \/3, the vector iTui is in the kernel of the stability 



matrix B so rhombus A becomes degenerate, which is in agreement with Theorem 2.7 since L = 



here. For rhombus B, we have y 2 < 3 — 2\/2 < 2 — for — 1 < m < 0, so the eigenvalues always 
form a real pair. 

Next we consider the eigenvalues associated to {t> 2 , Kvz}. In this case, the sign of u 2 — fj, 2 . is 
determined by the sign of the product g{y 2 )g{l/y 2 ). The discriminant of the cubic g is negative, so it 
has only one real root which we denote by k. Computing a Grobner basis for g(y 2 ) and equation ( 16 ) 
yields k = — m*(3m* + 2) « 0.1278. Moreover, substituting y 2 = k into equation (16) gives m = m*. 
The pair of eigenvalues corresponding to {v2,Kv2} are pure imaginary for k < y 2 < 1/k, real for 
< y 2 < k or y 2 > k, and when y 2 = k or 1/k. For rhombus A, we have n<l<y 2 <3 + 2y/2 < 
1/k, where the last inequality is rigorously shown by checking g(l/(3 + 2\/2)) > 0. Thus, this 
pair of eigenvalues is always pure imaginary for the rhombus A family. For rhombus B, we have 
/t<?/ 2 <l/3<l/fi;as long as — 1 < m < m*. For these m-values, the pair of eigenvalues is pure 
imaginary. At m = m* , y 2 = K and rhombus B becomes a degenerate relative equilibrium with v-i 
in the kernel of the stability matrix B. For m* < m < 0, y 2 < k and the eigenvalues bifurcate 
into a real pair. Combining these conclusions with those for the eigenvectors {v%, Kv{\ finishes the 
proof. □ 



Remark. 1. For completeness, we note that the equilibrium point for rhombus B at m = — 2+y/3 
is linearly unstable. This follows since the nontrivial eigenvalues of M~ 1 D 2 H(zq) are the two 
real pairs ±/ii and ±/x 2 - 



If m > 0, then Theorem |3.5| implies that the rhombus A family is actually nonlinearly sta- 
ble. This provides a possible explanation for the long-term persistence of the rhombus-like 
configuration in the numerical experiments of Kossin and Schubert [20J . 

For —2 + v^B < m < 0, rhombus A is linearly stable, but numerical calculations using Matlab 
show that it is a saddle of H restricted to I — Iq > 0. This shows that for opposite signed 
circulations, it is possible for a relative equilibrium to be linearly stable, but not a minimum 



of H restricted to a level surface of /. The proof of Theorem 3.2 breaks down because v{ Mv\ 



and v%Mv2 are each negative, so the value of q ZQ (v) can also be negative, leading to a saddle. 



4.3 The Isosceles Trapezoid Family 

Another one-parameter family of symmetric relative equilibria in the four- vortex problem is formed 
by isosceles trapezoid configurations. This family was described and analyzed in [14]. Set Ti = 
T 2 = 1 and T3 = T4 = m, treating m G 1R + as a parameter. Vortices 1 and 2 lie on one base of the 
trapezoid while vortices 3 and 4 lie on the other, and the two legs of the trapezoid are congruent. 
The pair of vortices with the larger circulations lie on the longer base, with m = 1 corresponding to 
a square. This family exists only for m > 0. As m — > 0, the configuration approaches an equilateral 
triangle with vortices 3 and 4 colliding in the limit. 

The coordinates of the trapezoid are given by £\ = (1,0), f 2 = (—1,0), £3 = (—x,y) and 
£4 = (x,y), which places the center of vorticity at c = (0,my/(m + 1)) (see Figure [3]). To be a 
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relative equilibrium, we must have 



(x, y) = (s/a, \Jlm + 3 — a) where a = 

These expressions were derived from the corresponding results in [14J. One can check that the 
positions Zi — fj — c satisfy the defining equations Q with angular velocity tu = (2m + 1)/ (2m + 2). 
This value of u concurs with the one obtained using the formula u = L/(2I). Note that x and y 
are both real only when m > 0. 
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Figure 3: The setup for the isosceles trapezoid family. Figure shown is for m = 0.3, with the center 
of vorticity indicated by the small solid disk. 



2m + 1 



Since the key matrix M^ 1 D 2 H(zq) is invariant under translation, we can evaluate it by substi- 
tuting in the values £{ given above. Due to symmetry, it is possible to give explicit formulas for the 
eigenvectors of M~ l D 2 H(z ). Note that rf 3 = r| 4 = (x + l) 2 + y 2 and r 2 4 = r| 3 = (1 — x) 2 + y 2 . 
The 2x2 off-diagonal blocks of D 2 H(zo) are given by 

-1 0" 
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y 2 -(x+ l) 2 -2y(x + 1) " 
-2y(x + l) {x + l) 2 -y 2 ' 
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y 2 -(l-a;) 2 -2y(l-x)~ 
-23/(1 -x) (l-x) 2 -j/ 2 ' 

Surprisingly, the null space of D 2 H(zq) is four-dimensional for all values of m. In addition 
to the two vectors coming from the conservation of the center of vorticity, the vectors wi,Kw\ 
are also in the kernel of D 2 H(zq), where W\ = [mx,0, —mx,0, 1,0, — 1,0] T . The vector w\ was 
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found by searching for vectors in the kernel of the form [a, 0, — a, 0, 6,0, — b, 0] , and was verified 



symbolically using Maple. Taking [i\ = in Lemma 2.4 shows that two of the nontrivial eigenvalues 
for the isosceles trapezoid family are Note that the invariant space formed by the span of 

w% and Kwi is M-orthogonal to the invariant space arising from the conservation of the center of 
vorticity. Thus, although the eigenvalues ±ui are repeated, the Jordan form of B does not have an 
off-diagonal block. 

The remaining nontrivial eigenvectors for M~ l D 2 H(z ) are harder to locate. After some initial 
numerical investigations using Maple, we found the pair of eigenvectors w 2 , Kw 2 with corresponding 
eigenvalues ±/z 2 , where w 2 = [—a, —m, a, —m, ax, 1, —ax, 1] 



r 



This, in turn, gives 



wiA/3(2m + 1) m 2 + m + 1 

— ^ and a 2 = —, r-, r. 

m 2 + m + l p 2(m + l)(m + 2) 



2 2 3(m 2 + 4m + 1) 3L 



4(m + 2) 2 4(m + 2) s 

It is interesting, though not altogether surprising, to see the reappearance of the total angular 
vortex momentum L here. It follows that the remaining nontrivial eigenvalues lie on the imaginary 
axis and the isosceles trapezoid family is stable for all m. We have proven the following result. 

Theorem 4.2. The isosceles trapezoid family is stable for all values of the parameter m > 0. The 

, ±V3L . 

nontrivial eigenvalues are ±ui ana — -i. 

y 2(m + 2) 



5 Conclusion 

We have adapted the approach of Moeckel from the n-body problem to study the linear stability of 
relative equilibria in the planar n- vortex problem. Due to some special properties of the logarithmic 
Hamiltonian in the vortex case, the stability is easier to determine and a complete factorization 
of the characteristic polynomial into quadratic factors exists when all the circulations have the 
same sign. Using a topological approach motivated by a conjecture of Moeckel's, we show that, in 
the case of same-signed circulations, a relative equilibrium is linearly stable if and only if it is a 
nondegenerate minimum of the Hamiltonian restricted to a level surface of the angular impulse. In 
this case, it follows quickly that linear stability implies nonlinear stability as well. Two symmetric 
examples of stable configurations in the four- vortex problem, a rhombus and an isosceles trapezoid, 
are analyzed in detail, with interesting bifurcations discussed. It is hoped that this work will be of 
some use to researchers in fluid dynamics, particularly to those who have discovered vortex crystals 
in their numerical simulations. 
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